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^sO , Using our previous approach to electromagnetic emission during tunneling, an ex- 

^ ! plicit, essentially classical, formula describing the bremsstrahlung spectrum in a decay 

| is derived. The role of tunneling motion in photon emission is discussed. The shape 

^ ■ of the spectrum is a universal function of the ratio E 7 /Eq , where E 7 is the photon 

energy and Eq is a characteristic energy depending only on the nuclear charge and the 
energy of the a particle. 



PACS numbers: 23.60.+t, 23.20.-g, 27.80,+w. 27.90.+b 



During the a decay of a nucleus the a particle tunnels through the Coulomb barrier 
^ ■ and is accelerated beyond the classical turning point to its final energy. Thus electro- 

magnetic radiation should be emitted during the process. In spite of the fundamental 
nature of the 7 emission duringa decay, the corresponding bremsstrahlung spectrum 
was neither measured, nor considered theoretically until recent years. Normally soft 
7 emission accompanying the Coulomb interaction of heavy particles is well described 
by classical electrodynamics. The case of a decay is somewhat special, since part of 
the 'trajectory' of the a particle lies within the underbarrier region and may be de- 
scribed classically only in some limited sense, as a motion in imaginary time. While 
the golden rule of quantum mechanics gives a straightforward recipe to calculate the 
emission probability, intuitively one would like to understand whether tunneling may 
be somehow incorporated in the general framework of classical electrodynamics and 
whether it makes any sense to say that photons are emitted during tunneling motion. 

In Ref. [1] the electromagnetic radiation by a charge tunneling through a potential 
barrier was considered under the conditions that the motion is quasiclassical, i.e. the 
barrier is smooth compared to the particle wave length, and that the energy of the 
emitted photons is small compared to the particle energy. It was shown that in this 
case the emission spectrum is described by the well known classical formula, involving 
the Fourier transform of the particle acceleration, with the only difference that the 
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Fourier integral should be taken along a special contour in the complex plane of the time 
variable. The resulting formula for the bremsstrahlung spectrum per tunneling particle 
is classical (does not contain Planck constant), even though it describes radiation 
emitted in the tunneling process, which is unknown to classical mechanics. The result 
of Ref. [1] may be directly applied to a decay. 

On the experimental side, D'Arrigo et al [2] were the first to report bremsstrahlung 
in a decay. Their data for 226 Ra and 214 Po gaveemission probabilities much larger than 
those predicted by the so-called Coulomb acceleration model, in which only acceleration 
outside the Coulomb barrier is taken into account, but close to those predicted by 
the sudden acceleration model. These results are strange since even the Coulomb 
acceleration model grossly overestimates the emission probability. Recently Kasagi 
et al [3] measured the 7 emission associated with the a decay of 210 Po, obtaining 
probabilities which are almost two orders of magnitude smaller than those of Ref. 2, 
and which have a quite different dependence on the photon energy. Kasagi et al applied 
the general formula of Ref. [1] and by numerical calculations obtained a good fit to 
their experimental data. However, numerical calculations for the specific case of 210 Po 
do not allow to understand the general features of the spectrum and its dependence 
on the nuclear charge and a particle energy. Thus, there is a need for an explicit 
theoretical formula describing the bremsstrahlung spectrum. 

Papenbrock and Bertsch [4] considered the problem and arrived to the conclusion 
that the main contribution to the photon emission stems from Coulomb acceleration, 
while the contribution of tunneling is negligible, contrary to the point of view expressed 
in Refs. [1] and [3] . Their numerical results also agree well with the experimental data 
of Ref. [3]. 

The purpose of the present Letter is to derive an explicit formula for the bremsstrahl- 
ung spectrum in a decay and to elucidate the role of tunneling. It will be shown that 
taking tunneling motion into account is very important to obtain the correct description 
of the bremsstrahlung. 

The a particle moves in the Coulomb potential U(r) = 2Ze 2 /r, which is truncated 
at the nuclear radius r . The classical turning point is defined by r c = 2Ze 2 /E a , where 
2e is the charge of the a particle, Ze is the charge of the daughter nucleus, E a = mv 2 /2 
is the a particle energy, m is the reduced mass, and v is the relative velocity at r = 00. 
We assume tq to be much smaller than r c . The initial state is described by an outgoing 
wave with the asymptotic form oc (l/r)exp(ifcjr), with fcj = (2mE a ) 1 ^ 2 /% . The final 
state at r = 00 has the form oc (l/r)cos(/c/r + 0), kf = [2m(E a — huo)} 1 / 2 /%, huo = E 1 
is the energy of the emitted photon, and is a phase, related to the scattering phase. 
In the underbarrier region the wave function of the initial state increases exponentially 
when r changes from r c inwards, while the wave function of the final state decreases 
from r c inwards. (For a detailed discussion of the initial and final states entering the 
matrix element for photon emission see Ref. [4]). 
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The following conditions will be assumed: 



hr c = 2Ze 2 /hv » 1, (1) 
fiuj « E a . (2) 

Eq. (1) says that the motion is quasiclassical which is generally very well justified, 
e.g. for Z = 50, E a = 4 MeV we have kir c = 100. Eq. (2) means that the energy losses 
are small, and we will see below that the main body of the bremsstrahlung spectrum 
lies in this energy range. 

We can now use the quasiclassical approximation for the initial and final wave 
functions by replacing kr with J k(r)dr , where ki = [2m(E a — [/(r)] 1 / 2 /h , kf = 
[2m(E a — fiw — [/(r))] 1 / 2 /h . The reflected part of the final state wave function with 
the asymptotic behavior oc exp(— ikfr) may be neglected [1], since its contribution to 
the matrix element is exponentially small, compared to the contribution of the outgoing 
wave oc exp(ikfr). To the first order in the small ratio hu/E a we have ki(r) — kf(r) = 
fuj/v(r), where v(r) is the classical velocity at point r. Proceeding further as in Ref. 
1, one finally obtains the classical formula [5] describing the bremsstrahlung spectrum: 



d£ _ 2 (Z eff e) 
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a w \ , (3) 



where d£ /dcu is the total energy emitted per decay per unit frequency range, Z e ff = 
2(A — 2Z)/(A + A) is the effective charge, A is the mass number of the daughter nucleus, 



/ a(t) exp(-iut)dt (4) 
J c 



is the Fourier transform of the acceleration a(t), the integral is taken along the contour 
C which is defined by the manner in which the time variable t = J r r c dr/v{r) changes 
when r runs from to oo along the real axis. In the underbarrier region t is imaginary, 
and it is easy to calculate the value of t for r = 0: t = —it, where r = (n/2)r c /v is the 
tunneling time. This value practically does not depend on the details of the potential 
curve at small distances and may be calculated by extending the Coulomb potential 
down to r = 0. In the following we will use the relation a w = iuv u , where is defined 
through the classical velocity v(t) analogous to Eq. (4). 

The function v w generally depends on the details of the potential function in the 
region close to the nuclear radius. However, these details are irrelevant for sufficiently 
low frequencies, such that ut\ « 1, where t\ is the characteristic time the particle 
spends in the vicinity of tq. This time can be estimated as t\ ~ tq/vq, with mv\ ~ 
Ze 2 /r , thus ti ~ T(r /r c ) 3 ^ 2 , and the condition ivti « 1 may be rewritten as cur « 
( r ' c / r o) 3 ^ 2 - Since, as shown below, the main part of the bremsstrahlung spectrum lies 
in the region cut ~ 1, and since r c /r >> 1, it can be seen that the structure of the 
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potential function in the vicinity of the nuclear radius may influence only the high- 
energy tail of the spectrum where the emission probability is already exponentially 
small. Thus, the main body of the spectrum may be calculated by extending the 
Coulomb potential down to r = 0. 

Introduce the dimensionless coordinate £ = r/r c , the dimensionless time n = vt/r c , 
and the dimensionless frequency v = (l/2)uor c /v = uoZe 2 /(E a v ). Then d£/doo may be 
expressed as 

/oo 
dr}(d£/dr})exp(—2ii/rf) (6) 
-ivr/2 

where the function £(77) is the solution of the (dimensionless) Newton equation d 2 £ s /drf = 
(2£ 2 )- 1 with the conditions £(0) = 1, (d£/dn) v=00 = 1. 

This may be compared with the results for two related classical problems: a) 
Coulomb acceleration model in which only acceleration outside the Coulomb barrier 
is taken into account. The lower limit for the integral in Eq. (6) should be n = 0, 
instead of n = —in/2, b) Head — on collision of two particles with charges 2e and 
Ze. In this case the lower limit in Eq. (6) should be n = —00, instead of n = —in/2. 

Introducing a new variable z, such that £ = (l+cosh,2)/2, n = (z+sinhz)/2, we 
obtain 

/oo 
dz sinh zexp[— %v{z + sinh 2)]. (7) 
-ITT 

The integration contour in Eq. (7) goes from —in to along the imaginary axis 
and from to 00 along the real axis. It may be safely deformed to go from —in to 
—in + 00 (see Ref. 5). Finally introducing x = z + in, we get 

roo 

J(v) = — iz/exp(— nv) / dxsinha:exp[ii/(smha; — x)\. (8) 

Jo 

By integration by parts J(v) may be also expressed in a form which is more suitable 
for numerical calculations: 

J(v) = exp(— 7rv)[l —iul(v)], 

where _ 

I{y) — / dx(l — e~ x ) exp[iz/(sinha; — x)\. (9) 
Jo 
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The integral in Eq. (9) converges more rapidly, than the one in Eq. (8). 

Eqs. (5) and (8) give an explicit formula for the bremsstrahlung spectrum in a 
decay, provided the conditions given by Eqs. (1) and (2) are fulfilled, and this is the 
main result of this work. For the head-on collision problem, mentioned above, J{v) 
should be replaced by Jcoiiiy), where J colli 1 ') ls gi ven by the same Eq. (8) but with the 
lower limit of integration replaced by — oo (see Ref. [5]). Thus J co u{y) = 2ReJ(v). 

The Coulomb acceleration model, which neglects tunneling, predicts a spectrum 
that falls down at high frequencies as uo~ 2 . The reason is that in this model the 
acceleration has a singular behavior, increasing step-wise from a zero value at t < to 
a finite value at t — 0. Taking tunneling into account makes the function a(t) smooth, 
and as a result, the spectrum falls down exponentially. One might say that there is 
a destructive interference between electromagnetic fields emitted during the tunneling 
motion at r < r c and during the classical acceleration stage at r > r c . (Similarly, 
in the classical Coulomb collision problem, the interference between the incoming and 
outgoing parts of the trajectory also produces an exponentially decreasing spectrum at 
high frequencies [5]). Thus, as can be also seen from Eqs. (4), (6), taking account of the 
tunneling motion is very important to obtain the correct description of bremsstrahlung. 

The opposite conclusion reached by Papenbrock and Bertsch [4] is due to a mis- 
interpretation. When applied to the quasiclassical limit of the acceleration matrix 
element (Eq. (8) above), their reasoning goes as follows. Noticing that lmJ(u) = 
at v — [6], they make the statement that ImJ(z/) is numerically small compared to 
ReJ(u) for arbitrary v (in fact, for v » 1, the ratio Im J(u)/Re J(u) is equal to — 1/\/3 
). On these grounds the imaginary part of J{v) is discarded, while the real part, as 
shown above, can be expressed as (1/2) J co u[y). Since calculation of J co u{y) does not 
involve tunneling motion, the conclusion is made, that the contribution of tunneling 
is negligible. In fact, neglecting ImJ(z/) is equivalent to replacing the true initial wave 
function oc exp(i/cjr) at r = oo) by its real part which describes an a particle scattered 
by the nucleus (oc cos(kir) at r = oo), but normalized to 1/2 incoming flux. Not 
surprisingly, the spectrum for a decay, obtained in Ref. 4 after neglecting lmJ(p), is 
equal to 1/4 of the spectrum for the head-on Coulomb collision problem [7]. However, 
this result has a correct behavior both at low and at high frequencies, the difference 
with the true spectrum being not very large (25% for v » 1), especially when plotted 
on a logarithmic scale. Clearly, the fact that the spectrum for the collision problem 
(with a factor 1/4) does not differ very much from the true spectrum in a decay, does 
not mean that tunneling is not important. 

Papenbrock and Bertsch [4] have also performed a numerical calculation of the 
exact matrix element, without neglecting its real part and without using the quasi- 
classical approximation, for the specific case of 210 Po. The resulting spectrum is very 
close to the spectrum given by the explicit Eqs. (5), (8), (9) of the present work (see 
below). Thus the numerical results of Ref. 4 are correct and are in agreement with 
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the above formulae, which explicitly take into account the tunneling part of the tra- 
jectory. However the physical interpretation of these results in Ref. 4, leading to the 
conclusion that the main contribution to the photon emission stems from the Coulomb 
acceleration, is erroneous. In fact, it is the contribution of the tunneling region that 
accounts for the very large difference between the true spectrum and the predictions 
of the Coulomb acceleration model. 

It is easy to find the asymptotic form of J(v) both for small and for large values of 
v . For v « 1 J(v) = l-7rz//2-?z/ln(27/V), where 7 =1.781, so that \J(v)\ 2 « l-nu, 
while for v » 1 



J(y) = -iexp(i7r/3)r(2/3)(4z//3) 1/3 exp(-7rz/) (10) 



\J( V )\ 2 = r 2 (2/3)(4z//3) 2/3 exp(-27n/) = 2.22z/ 2/3 exp(-27rz/). (11) 

It follows from Eqs. (10), (11) that for v » 1 | J coU {v)\ 2 = 3| J{v)\ 2 . Fig. 1 shows the 
normalized spectrum, given by the function |J(z/)| 2 , calculated numerically, together 
with the asymptotic curves for v « 1 and v » 1. It can be seen that Eq. (11) 
becomes a very good approximation for v > 0.5. In this region the |J(^)| 2 depen- 
dence is dominated by the exponential factor in Eq. (11), which may be expressed as 
exp(— E 7 /E )), where 

is a characteristic energy, 2nu = E^/Eq. For the experimental conditions of Ref. [3] 
(Z = 82, E = 5.3 MeV) one finds E m 75 keV. Note that the shape of the spectrum 
is a universal function of E 7 /Eq. 

The experimentally measured quantity is the differential probability of photon emis- 
sion per decay d 2 N/dE 1 d£l = (AnhE^dS /du. Thus 



d 2 N _ Zl f} e 2 E a 1 
dE^dVt 37r 2 Tic mc 2 E^ 



\J{u)\\ (12) 



This quantity is plotted in Fig. 2 as a function of E y for the case of a decay of u Po 
(heavy line), together with the experimental data taken from Ref. [3]. The result of the 
numerical calculation in Ref. [3] is presented by the dotted line. This calculation was 
done by using the general formula of Ref. [1] , and thus should give the same result as 
the present work. However, one can see that the two results are qualitatively different: 
a)the slopes at low 7 ray energies differ by a factor ~ 2 and b)while the present theory 
gives a spectrum that falls down exponentially at high E 1 , the calculation of Ref. [3] 
gives a broad maximum around E~ ~ 500 keV. The difference in the slopes can not be 
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easily understood. As to the high energy part of the spectrum, we should note that 
there is a limitation of the present theory related to our neglecting the details of the 
inner barrier region and the extension of the Coulomb potential down to r = 0. As 
explained above, this is justified if E 1 « £ (r c /r ) 3 / 2 860 keV (r c /r ps 5 for 210 Po). 
For energies comparable to or higher than this value the shape of the spectrum should 
depend on the parameters of the inner barrier region. Thus, the maximum on the curve 
calculated in Ref. [3] (where the cutoff of the Coulomb potential at r = r was taken 
into account) might be related to this circumstance. More precise experimental data 
are needed to determine unambiguously the true slope at low energies, as well as the 
existence of any specific features at the high energy tail of the spectrum. 

In summary, an explicit, essentially classical, formula, describing the bremsstrahlung 
spectrum in a decay, was derived, and it was shown that taking the tunneling motion 
into account is crucial for the correct description of the spectrum. The shape of the 
spectrum is a universal function of the ratio E 1 / E Q , where E 1 is the photon energy and 
a Eq is a characteristic energy depending only on the nuclear charge and the energy of 
the a particle. 

I thank Michael Shur for his hospitality during my stay at Rensselaer Polytechnic 
Institute, where this work was completed. I appreciate useful discussions with Thomas 
Papenbrock. 
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Figure captions 

Fig. 1. Normalized bremsstrahlung spectrum, (a) - normal scale: 1 - numerical 
calculation, 2- low frequency asymptote |J(^)| 2 = 1 — iru, 3 - high frequency asymp- 
tote, Eq. (11). (b) - logarithmic scale: 1 - numerical calculation, 3 - high frequency 
asymptote, Eq. (11), 4 ~ result for Coulomb acceleration model, 5 - high frequency 
asymptote for Coulomb acceleration model |J(^)| 2 = (4z/) -2 . 



Fig. 2. 7 ray emission probability in a decay of 210 Po. Results of present theory 
(heavy solid line) are presented together with the experimental data points and results 
of numerical calculation from Ref. [3] . Thin solid line is the high frequency asymptote 
given by Eqs. (11), (12). 
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